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I use the Newtonian equation of hydrostatic equilibrium for an isotropic fiuid sphere to generate 
exact anisotropic solutions of Einstein's equations. The input function is simply the density. An 
infinite number of regular solutions are constructed, some of which satisfy all the standard energy 
conditions. Two classes of these solutions generalize the Newtonian polytropes of index and 1. 
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I. INTRODUCTION 

The study of exact spherically symmetric perfect fluid solutions of Einstein's equations has been reinvigorated 
over the last decade or so and we now have a remarkable collection of solution generating techniques At the 
same time, interest in exact anisotropic fluid solutions of Einstein's equations has grown [2]. Whereas a solution 
generating technique is known for anisotropic fluids [3] , this requires the specification of two input functions and one 
is not guaranteed a reasonable output, say in the density p. A less general technique is given here, but the algorithm 
presented requires only one input function, the density itself. The technique is demonstrated by the construction of 
an infinite number of regular solutions, some of which satisfy all the standard energy conditions. 

II. GEOMETRY 

Let us write the geometry in the familiar form 

^ 5 + r^dn" - e2*('-)dt2 (1) 

1 — 2m[rj/r 

where dH.^ is the metric of a unit 2-sphcrc {dd'^ + sin^9d(p'^). The source of (H]), by way of Einstein's equations, is taken 
to be a comoving fluid described by the stress-energy tensor — diag[p{r), P{r), P{r), —p{r)]. Einstein's equations 
give the effective gravitational mass 



m = An / p{x)x'^dx, (2) 







the source equation for $ 



$' ^ + (3) 
r{r - 2m) ' ^ ' 



and the generalized Tolman-Oppenheimer-Volkoff equation which we write in the form 

P=^{p' + ip + p)^')+P, (4) 

where ' = d/dr. The fluid distribution is assumed to terminate at r = i? where p{R) = and to join there, by way of 
a regular boundary surface, onto a Schwarzschild vacuum of mass m(R) = M. No further conditions need be specifled 
for this junction. 



A fluid is considered regular if it is free of singularities. Scalar polynomial singularities involve scalar invariants built 
out of the Riemann tensor. Due to the spherical symmetry assumed here, there are only four independent invariants 
of this type [sj . A direct calculation shows that p and P do not enter these invariants and that the Ricci invariants 
are regular as long as p,p and P are. However, the second Weyl invariant grows like 

rp +2{p-P) \ !_ 
p-3p + 4:P J ^ ' 
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as r — > 0. As a consequence, we have the following necessary conditions for regularity 

p(0)-0, p{0)^PiO), (6) 

the latter of which is already obvious from One could, of course, incorporate derivatives of the Riemann tensor 
into the construction of invariants. For example, a calculation shows that DTZ grows like 

^ (7) 



as r ^ where TZ is the Ricci scalar and □ is the covariant d'Alembertian. However, the background theory involves 
the Riemann tensor, not its derivatives, and so we impose no a priori restrictions beyond ([6]). 



III. THE ALGORITHM 



The algorithm proposed here involves the use of the Newtonian hydrostatic equilibrium equation for an isotropic 
fluid, 

P = ^ — - (8) 

In effect then p and p are related exactly as in the Newtonian theory of isotropic fluids, e.g. [6|. A solution in 
Einstein's theory is obtained simply by reading off P from (|4]). No isotropic solutions for Einstein's theory can be 
found in this way since if we set P = p in (jl]), and use ([8]), the resultant quadratic equation for p returns p{0) < 0. 
One can take the view that it is the tangential stress P that holds these configurations together. From dU and ([5]), 
and assuming p > 0, we have 

P > p > 0. (9) 

Explicitly, inserting ([3]) and jS]) into ([U we obtain 

Imp 1 (p + p){m + Anr^p) 



2 r 2 r~2m 



(10) 



As a result, the weak and classical strong energy conditions are automatically satisfied via the proposed algorithm, as 
long as p > 0. The dominant energy condition is not ?|. However, in the present circumstance, the classical strong 
energy conditions are in fact weaker than the weak energy conditions and so the trace energy condition 

2P + p<p (11) 

is usually applied. Subject to (fTTj) it is known that the Buchdahl bound R/2M > 9/8 f§\ is robust ^. From (dll) we 
can write 

P , . 

- < e (12) 

and consider the dominant energy condition as e = 1 and the stronger trace energy condition as e = 1/2 [loj . 



IV. EXAMPLES 



The procedure outlined above can be executed with the specification of p alone. We demonstrate this here with an 
infinite number of simple examples mostly derived from the ansatz 

p = -ar" + 6 (13) 

where a and b are constants > and n is an integer > 1. First, let us consider the simplest case, that of constant 
density, a — [ill] . 



3 



A. a = 

From ^ and ^ it follows that m = Airbr^/S and so M = ATrhR^/i. From (Hj), and setting p{R) = 0, we obtain 

p=\TTh^ {R-r){R + r) (14) 

and so p = — 47r6^r/3. This is a Newtonian polytrope of zero index, e.g. We can now simply read off P from ^ 
with dg. We find 

P= , "^^^^ (r2(47r26V-3 + 8^6r2)-i?2(8^2^2^4„4^2^2^2^2_3)N ^^g) 
3(3 - 87r6r2) ^ ^ ' ^ " ^ ' 

To avoid singularities in P then 



and so for non-singular solutions R < 2M as expected. At r = we find 

p(0) = 6, p(0)= p(0) = P'(0) = (17) 



along with 



At r = i? we find 



along with 



„0) = P(0) = H^?5!. (18, 



p{R) ^h, p {R) = p{R) = 0, p (i?) = — i? (19) 



(20) 



For solutions with P{R) < p, since P{r) > 0, the dominant energy condition is also satisfied throughout the solution. 
From pO)) for these solutions we find 

where P{Rc) — p, and so in the limiting case Rc/2M — 1/2(^/2— 1) ~ 1.207 somewhat larger that the the Buchdahl 
bound R/2M = 9/8 = 1.125 iJ]. Some examples are shown in Figure [T] In order to demonstrate that energy 
conditions of the type (I12p can be satisfied, it is convenient to define = br^ and = bR^ so that 



^ ^'^ (x2(4 7rV- 3 + 871x2) -X2(87r2x^-47rVx2- 3)). (22) 



p 3(3-87rx2) 
A number of examples are shown in Figure [2l 



o > 



We set 



a 



l/n 



(23) 
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FIG. 1: The variables P (increasing) and p (decreasing to at r = i?) for the critical case R — Rc and for R = Rc — 0.05 (for 
which R/2M ~ 1.707) for the constant density solutions with 6 = 1. 



SO that p{R) = 0. Other boundary conditions are given below. From ^ and ([T^ it follows that 

_ 47r 3 /(3 + n)6- 3ar" 



and so 



From dHl), ini) and ^ we obtain 



47r „o / nb 
M = — i?3 

3 \3 + n 



(24) 



(25) 



_ 27r / -r2 (62 (3 + n) (2 + n) {1 + n) + ar"- (3 ar" (2 + n) - 2 6 (6 + n) (1 + n))) + R%'^n^ (4 + n) \ 
^"TI (3 + n) (2 + n) (1 + n) j' ^^'^^ 



We can now simply read off P from ^ with ([3]). We find 

P= 27r \ ( H{r) 

\2,{i + n){2 + nf{l + nY) \V{r) 

where is a polynomial of degree 6 + 4n, given explicitly in Appendix A, and 

■P = 3(3 + n) - 87r(3 + n)hr'^ + 247rar2+". 

As shown in Appendix B, (I28|) has no roots for r > as long as 

'3 + n\ / 167r6 \"/^+V^\"/2 

a > ■ 



Stt y V3(2 + ?i) 



(27) 



(28) 



(29) 



Subject to condition then the solutions given here are free of singularities throughout the range < r < i?. 
Whereas the configurations satisfy the weak and strong energy conditions, the dominant energy condition is satisfied 
for sufficiently large a. This is demonstrated below. 




FIG. 2: The ratio P/ p for models characterized hj X = Xc — S where 5 = for the top curve and 5 increases by 0.01 for each 
subsequent curve. The insert shows more detail. For these models R/2M = 3/87r(Xc ~ (5)^. 



At r = wc find 



p(0) = 6, p'(0) = -anr""\ p'(0) = F'(0) = 



(30) 



along with 



3 ) V(3 + n)(2 + n)(l + n) 



(31) 



Ax, r = R we find 



p{R) - p{R) = P{R) ^piR)^0, p (R) 



nb 



(32) 



along with 



P{R) 



16 / R^n%^n^ 

y V (-3(3 + n) + 8 TT 6i?2n)(3 + n) 



(33) 



From (USD, ^ and dMl) we find 



R _3(3 + .)_^ /3W2(3+^y/"-^^^ ^3^^ 



2M 87m&i?2 V3(2 
Since lim a — I, the solutions given here allow an arbitrarily large degree of compactness, in principle. However, 

n — ^oo 

smaller values of a correspond to smaller values of a which lead to a violation of the dominant (and therefor trace) 
energy condition. Again, this is demonstrated below. From it follows that the Buchdahl bound R/2M = 9/8 is 
obtained for 

/37rn62/"+i\"/2 



For n — 1 there is no "core" in the density. It follows from ((29)) that 

a > ^V2^6^/^ (36) 

and from (j34[) that R/2M > 256/243. The physical variables are shown in Figure [3] for a = 2/3\/2^63/2 and in Fi eure 
IHfor a — ■y/3/8\/27r&'^/^, the Buchdahl limit. In the first case all energy conditions (excepting the trace) are satisfied 
and a = 4/3. (The trace energy condition can be satisfied with a slightly larger a, for example, a = I9/27-\/27r6'^/^.) 
In the second case the dominant energy condition fails. Unlike the case a = 0, the minimum a for which the dominant 
energy condition is satisfied must be found numerically for each n. For n = 1 the minimum a ~ 1.29, of course above 
the Buchdahl limit. 
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FIG. 3: The variables p,P,p and m for a = 2/3\/2Tvb'^^^ and b = 1 for the linear density solution. 



2. n = 2 

For 71 = 2 it follows from that 

a > ^TTb', (37) 

and from that R/2M > 25/24. The physical variables are shown in Figure [S] for a = 13/97r6^ and in Figure 
[6] for a = 6/57r6^, the Buchdahl limit. In the first case all energy conditions (excepting the trace) are satisfied and 
a — 65/48 ~ 1.354. (The trace energy condition can be satisfied with a slightly larger a, for example, a — 14/97r6^.) 
In the second case the dominant energy condition fails. The solutions with n > 2 are qualitatively similar, but with 
a flatter density profile for smaller r, that is, a larger "core". An example is shown in Figure [7l 

3. Polytrope of index 1 

The Newtonian polytrope of zero index was discussed above, but it is well known that there are two finite analytical 
polytropic solutions, e.g. [61]. In this section we use a classical polytrope of index 1 to generate an anisotropic fluid. 
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FIG. 4: The variables p, P, p and m in the Buchdahl limit a = •^3/8%/27r6'''''^ and 6 = 1 for the linear density solution. 
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FIG. 5: The variables p, P, p and m for a = 13/97r6^ and 6 = 1 for the quadratic density solution. 



Let us take 



P = Po- 



sin(ar) 



(38) 



where a is a constant, so that from ([5]) 



27r 2 
V= -tP 



(39) 
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FIG. 6: The variables p, P,p and m in the Buchdahl hmit a = 6/57r6^ and b — 1 for the quadratic density solution. 
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FIG. 7: The variables p,P,p and m for a ~ SSyr'^fo^ and 6 = 1 for the solution with n = 16. This case gives the minimum 
a ~ 1.24 for which the dominant energy condition is satisfied. 

and 

R^-. 40) 
a 

From ^ and ([38]) we have 

4:77 PQ /,s 

m = — 5— (sin(ar) — ar cos(ar)) (41) 
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and so 



from which wc have 



47r2 4 , 



R a2 



2M Sirpo 

and so > 9TTpg for the Buchdahl bound to hold. We now find 



P = 



I{r) 



32 7rr5a8 (r - 2 m) 

where / is a polynomial of degree 6, given explicitly in Appendix C. At r = we find 

p(0) = po,p(0) = P(0) = %pI p (0) = p (0) - P'(0) = 0. 



(42) 



(43) 



(44) 



(45) 



At r = i? we find 



along with 



p{R)^p{R)^P{R)^p{R)^0, p\r) = -P^ 



P{R) 



Wirpl 



a{a? — Sttpo) 



< 



(46) 



(47) 



where the last inequality assumes the Buchdahl bound. These models are qualitatively similar to the case n — 2 
discussed above. In a manner analogous to Figure [5] some details associated with these models are shown in Figure 

m 




FIG. 8: The ratio P/ p for models characterized by a = (3 + 5)yjTipo where the top curve has 5 = (the Buchdahl bound). 
Going down the curves, 5 increases by 0.1 except for the fourth curve down for which 5 = 0.25 and the bottom 3 curves for 
which 5 = 1.6, 2 and 2.5 respectively. For smaller values of 5 the curves are truncated. The insert shows more detail. 
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V. 



DISCUSSION 



An algorithm which converts isotropic Newtonian static fluid spheres into anisotropic Einsteinian static fluid spheres 
has been given which requires one input function, the density. An infinite number of explicit regular solutions have 
been generated, some of which satisfy all the standard energy conditions. 
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H{r) = b^n^R^in + 4)[4:R^b'^r'^n'^n^{n + 4) - 8 b^r*7r^{n + 3)(n + 2){n + 1) + 16 6r"+Va(n + 6)(n + 1) (Al) 

~6 7i7rr"+2a(n + 2){n + 1) - 247rV+2"a2(7i + 2) + 3 (ti + 3)(n + 2){n + 1)] 
+b^r^{n + 3)2 (n + 2f{n + 1)^(4 ^Vtt^ -3 + 8 br'^n) + 6 abr''+^{n + 3)(n + 2)(n + 6)(7i + 1)^ 
-2 7r62ar"+4(n + 36)(?i + 3)(n + 2f{n + if - 16 n^b^ar^'+^in + 3)(n + 2)(n + 6)(7i + if 
-9a2r2+2"(n + 3)(n+ l)(n + 2f + 12 a^nbr'^+^'^n + 2){3n^ + 20n + 36)(n + 1)^ 
+8n^b^a^r^+^''{n+ l){5n^ + 47n^ + 144n+ 108) - 18 aVr''+^ "(4 + 3 n)(n + l)(n + 2)^ 

-48 n\^br^+^'\n + 6)(n + 2)(n + 1) + 36nW^+^'\n + 2f . 



Acknowledgments 



possible by use of GRTensorll p^ . 



APPENDIX A: H 



For a > and n > 1 we find 



APPENDIX B: PROOF OF ([29]) 



This follows immediately from the following elementary theorem [12]: For 



1 - ax"^ + bx 



(Bl) 




(B2) 



APPENDIX C: / 



In terms of §Q and JHI) we find 



/(r) = 256^V {po~pf{po + pf 
128 TT^p r^a^ (aV^ {3 p^ + po^) - 4r {po - p) {po + p) + 2m {po - p) {po + p)) 
32 TrV^a** (6 a^mp^r^ + 2 p^r^ - 2 p^fr^ - 5 pfmr + 13 mp^r - 3 p^ + p^frn^) 

+8 7rp rara^ (3 c?mr' + 4 + 14 mr - 2 m^) 
—a m (4r +10TOr — m ). 
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